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Resumen 


Saucedo, H., Zavala, M., €z Fuentes, C. (septiembre-octubre, 
2015). Diseño de riego por melgas empleando las ecuaciones 
de Saint-Venant, y Green y Ampt. Tecnología y Ciencias del 
Agua, 6(5), 103-112. 


Se presenta un método de diseño de riego por melgas basado 
en la aplicación de un modelo que emplea las ecuaciones de 
Saint-Venant para describir el flujo del agua sobre el suelo, 
y la ecuación de Green y Ampt para representar el flujo del 
agua en el suelo. La principal ventaja del modelo presentado 
es el poco tiempo de cómputo que requiere para su 
aplicación, en comparación con el necesario para ejecutar un 
modelo hidrodinámico completo reportado en la literatura, 
mismo que emplea las ecuaciones de Saint-Venant para el 
flujo del agua sobre el suelo, acopladas internamente con la 
ecuación de Richards, que permite modelar el flujo del agua 
en el suelo. Con base en los resultados obtenidos mediante 
la aplicación del modelo para diez tipos de suelo diferentes 
se han elaborado cuadros de diseño operativo del riego por 
melgas. 


Palabras clave: riego por melgas, ecuaciones de Saint- 
Venant, ecuación de Green y Ampt, gasto óptimo. 


Introducción 


El objetivo del diseño de riego por melgas es 
conseguir la aplicación de una lámina de riego 
requerida por el cultivo de la manera más 
uniforme posible, conservando una eficiencia 
de aplicación alta. El diseño del riego consiste 
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en determinar el gasto de aporte y el tiempo 
durante el cual se aplica dicho gasto en la 
cabecera de la melga para lograr la mayor 
uniformidad posible; es decir, en determinar 
el gasto óptimo para una longitud de melga y 
una caracterización hidrodinámica del suelo 
específica. 
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Rendón, Fuentes y Magaña (1997) obtuvie- 
ron que el gasto óptimo de diseño guarda una 
proporción lineal con la longitud de la melga en 
la cual debe aplicarse. El resultado se obtiene 
usando un modelo formado por las ecuacio- 
nes de Lewis y Milne (1938) para describir el 
flujo del agua sobre la superficie del suelo, y 
de Green y Ampt (1911) para describir el flujo 
del agua en el suelo. Sin embargo, debe tenerse 
en cuenta que las ecuaciones de Lewis y Milne 
(1938) son la forma mas simplificada de modelar 
el flujo del agua a superficie libre en el riego por 
melgas. 

El objetivo del presente estudio fue desa- 
rrollar un modelo donde se utilice una forma 
precisa de modelar el flujo del agua a superficie 
libre, para lo cual se hace uso de las ecuaciones 
de Saint-Venant, mientras que para la descrip- 
ción del flujo del agua en el suelo se utiliza la 
ecuación de Green y Ampt, que es una forma 
sencilla de modelar el fenómeno manteniendo 
bases físico-matemáticas en la representación. 


Flujo del agua sobre la superficie del 
suelo 


El flujo del agua con una superficie libre es 
modelado con las ecuaciones de Saint-Venant, 
que resultan de la aplicación de las leyes de con- 
servación de masa y cantidad de movimiento. 
En una melga, la relación entre su ancho y el 
tirante de agua permite considerar las ecuacio- 
nes correspondientes al escurrimiento sobre una 
superficie de ancho infinito (Woolhiser, 1975). La 
ecuación de continuidad se escribe como: 


To (1) 
ot 0x ot 


la ecuación de momentum se escribe en la forma 
recomendada por Saucedo, Fuentes y Zavala 
(2005): 


197 2q 0 q oh 
+ +|2- 
hot h ox Max 


_ qa 2 
+2(] 1)+B er 0 (2) 
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donde q(x,t) = Ul(x,t)h(x,t) es el gasto por 
unidad de ancho de melga [L?T"]; x, la coor- 
denada espacial en la dirección principal del 
movimiento del agua en la melga [L]; t, el 
tiempo [T]; U, la velocidad media; h, el tirante 
de agua [L]; J, la pendiente topográfica de la 
melga [LL"]; J, la pendiente de fricción [LL”]; 
V, = 01/ ot, el flujo de infiltración [LT"], es decir, 
el volumen de agua infiltrado en la unidad 
de tiempo por unidad de ancho y por unidad 
de longitud de la melga; I, la lámina infiltra- 
da [Ll]; g, la aceleración gravitacional [LT?]; 
p = U, /U, donde U, es la proyección en la direc- 
ción del movimiento de la velocidad de salida 
de la masa de agua debido a la infiltración. 

La relación entre las variables hidráulicas q 
y h con la pendiente de fricción, denominada 
ley de resistencia hidráulica, es adoptada de 
acuerdo con Fuentes, De León, Saucedo, Par- 
lange y Antonino (2004), haciendo uso de una 
ley potencial de resistencia: 


ao ] (3) 
v 


donde v es el coeficiente de viscosidad cinemáti- 


ca del agua [L?T"] y k es un factor adimensional; 
d es un parámetro adimensional, que varía entre 
1/2<d<1, en función del tipo de flujo; los va- 
lores extremos d = 1/2 y d = 1 corresponden, 
respectivamente, al régimen de Chézy y al flujo 
laminar de Poiseuille. 

Para establecer las condiciones inicial y de 
frontera que deben sujetar a las ecuaciones de 
Saint-Venant para modelar el riego por melgas 
se tiene en cuenta una melga cerrada, para 
evitar la pérdida de masa fuera del dominio de 
riego. Para la fase de avance se tiene: 


g(x,0)=0 y  h(x,0)=0 (4) 


q(0,t)=4,, g(x,,t)=0 y h(x,,t)=0 (5) 


donde x(t) es la posición del frente de onda para 
el tiempo t y q, es el gasto de aporte en la entra- 
da de la melga. Para la fase de almacenamiento: 
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q(0,t) =q,, q(L,t)=0 (6) 
donde L es la longitud de la melga y q, es el 
gasto de aporte en la cabecera de la melga. Para 
la fase de consumo: 
q(0,t)=0, g(L,t)=0 (7) 
Para cerrar el sistema es necesario conocer la 
forma en que evoluciona en el tiempo la lámina 
infiltrada en toda posición sobre la melga, esto 
es, la ley de infiltración. 


Flujo del agua en el suelo 


El modelo de Green y Ampt (1911) se establece 
a partir de la ecuación de continuidad y la ley 
de Darcy con las siguientes hipótesis: a) el perfil 
de humedad inicial en una columna de suelo 
es uniforme O = 6 ; b) la presión del agua en la 
superficie del suelo es hidrostática: y =h > 0, 
donde ) es el tirante de agua; c) existe un frente 
de humedecimiento bien definido, caracterizado 
por una presión negativa: y =-—h,< 0, donde h, 
es la succión en el frente de humedecimiento; 
y d) la región entre la superficie del suelo y el 
frente de humedecimiento está completamente 
saturada (flujo en pistón), 8 = 0, y K=K_,, donde 
K_ es la conductividad hidráulica a saturación, 
es decir, el valor de la conductividad hidráulica 
de la ley de Darcy correspondiente al contenido 
volumétrico de agua a saturación. La ecuación 
diferencial ordinaria resultante es la siguiente: 


(h,+h)A0 
+ A 


¿1 (8) 


donde A6 = 6. - 6, es la capacidad de almacena- 
miento; Í es el volumen infiltrado acumulado 
por unidad de superficie de suelo o lámina 
infiltrada. 


Solución numérica 


Para el cálculo del gasto óptimo se hace uso de 
un modelo numérico para las cuatro fases del 


riego por melgas desarrollado por Saucedo, 
Zavala y Fuentes (2011). La disposición de las 
celdas de cálculo para las fases de avance, al- 
macenamiento, consumo y recesión se muestra 
en la figura 1. 

La forma discreta de la ecuación de conti- 
nuidad para la fase de avance se escribe como: 


) 


] 


[oq, +(1-0)g, ]ót-(x, -x 


o(h, +1,) 
+(1-0)(%, +1,) 


h +1 
oy oo. Jurt, died 


+(1-0)5, +1,,) 


-[6h, +(1-9)h, +01, +(1-0)1, |(x, -x,) 
+[9h, +(1-9)h, +91, +(1-9)1, ](x, -x,)=0 (9) 


para las fases de almacenamiento, consumo y 
recesión se tiene: 


n 


Figura 1. Disposición de las celdas de cálculo para las 


fases de avance, almacenamiento, consumo y recesión. 
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[v7, +(1- 0)q, 5 -[03, +(1-0)g,, 9 
-[9h, +(1- Q)h, +91, +(1- 9)1, ](x, -x,) 


+97, +(1- 9), +01, +(14)1,, Gx, -x;)=0 (10) 


La ecuación de momentum guarda la misma 
forma discreta para las cuatro fases del riego: 


27h(g, -q,)0t+(gh? 7), -h,)0t 


+h(x, -x,)[09, +(1 0), -wq,, -(1 0)q,] (11) 


+99h*(x, -x,)[0], +(1-0)J,-J,] 


+porgh(x, -x,) ol, +(1-0)1, -ol,, - 1-0), |=0 


En las ecuaciones (9) a (11), Ót es el paso de 
tiempo, w y q son factores de peso en espacio 
y tiempo, y las demás variables se definen 
de acuerdo con la figura 1. El cálculo de los 
coeficientes se realiza con base en los valores 
pertenecientes al nivel de tiempo anterior: 
q =(1-),, + 69, h=(1- 9), + 4h, (Saucedo et 
al., 2005). Se introducen las variables pequeñas 
(Strelkoff €: Katopodes, 1977): 0h, 5q,, 0h, y 9q,, 
tales que, para las celdas interiores se tiene: 
h, = h, + 0h, h,=h,, + 6h, q, = q,+ 99, y q,= 9, 
+ 0q,. Se sustituyen en las formas discretas de 


las ecuaciones de continuidad y momentum, y 
se obtiene un sistema de ecuaciones algebraicas, 
cuya solución permite avanzar en el tiempo los 
valores de las variables hidráulicas. El detalle 
de las entradas al sistema matricial puede con- 
sultarse en Saucedo et al. (2011). 

La ecuación de Green y Ampt (8) se resuelve 
numéricamente usando un método de diferen- 
cias finitas. El procedimiento se encuentra bien 
documentado en la literatura y puede consul- 
tarse, por ejemplo, en Burden y Faires (1985). 

Se ha utilizado un paso de tiempo constante 
At =1.0 s para el acoplamiento de las ecuaciones 
de Saint-Venant, y Green y Ampt. La discretiza- 
ción utilizada para la solución de la forma com- 
pleta de las ecuaciones de Saint-Venant guarda 
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semejanza con las reportadas en la literatura: 
Katopodes y Strelkoff (1977): At, =5 s; Akanbi 
y Katopodes (1988): At_. =1 s; Playán, Walker 


máx 


y Merkley (1994) At_. =2.12 s. 
Ventajas del modelo simplificado 


La principal ventaja que se tiene con el empleo 
del modelo simplificado basado en el uso de 
las ecuaciones de Saint-Venant, y Green y 
Ampt es el menor tiempo de cómputo que se 
requiere para la obtención de resultados de la 
aplicación del modelo. A manera de ejemplo, se 
utilizan datos de un experimento reportado en 
la literatura por Fuentes (1992). Se dispone de 
los siguientes datos: gasto unitario q, = 0.0032 
m*/s/ my; pendiente topográfica J, = 0.002 m/m; 
longitud de la melga L= 100 m,; parámetros para 
la ley de resistencia de Fuentes et al. (2004): 
d =1, k= 1/54; parámetro en la ecuación de 
cantidad de movimiento: fB = 2; valor inicial 
del contenido volumétrico de agua 0, = 0.2749 
cm*/cr”. Los parámetros 0, 0, m y n de las 
características hidrodinámicas fueron deter- 
minados a partir de la curva granulométrica 
del suelo franco de Montecillo. Siguiendo la 
metodología sugerida por Fuentes (1992) resulta 
9. = 0.4865 cm*/cm?, 8, = 0.0 cm*/cm*, m = 1258, 
n=2.2878 y n=11.0. Los valores delos parámetros 
y, = 32.75 cm y K_=1.84 cm/h fueron obtenidos 
mediante la aplicación de un procedimiento 
inverso para describir datos de una prueba de 
riego efectuada en el experimento mencionado. 

En una computadora de escritorio, con 
procesador Intel Y Core *M 17-4770 CPU O 3.40 
GHz y 8.0 GB de memoria RAM, la aplicación 
del modelo hidrodinámico completo basado 
en el uso de las ecuaciones de Saint-Venant y 
Richards para modelar las cuatro fases del riego 
por gravedad (avance, almacenamiento, consu- 
mo, recesión) requiere de un tiempo de cóm- 
puto de 3 505 segundos (58 minutos y medio, 
aproximadamente), mientras que la modelación 
empleando el modelo simplificado basado en el 
uso de las ecuaciones de Saint-Venant, y Green 
y Ampt requiere de sólo 23.6 segundos; es decir, 
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en el ejemplo bajo análisis, el modelo simpli- 
ficado es 150 veces más rápido que el modelo 
hidrodinámico completo. La comparación de los 
resultados obtenidos durante la simulación de 
la fase de avance empleando ambos modelos 
se muestran en la figura 2 para el caso de la 
posición del frente de avance, y en la figura 3 
para el caso de la evolución en el tiempo del 
tirante y la lámina infiltrada en la cabecera de 
la melga. Puede apreciarse que, no obstante la 
muy significativa diferencia en los tiempos de 
cómputo necesarios para la aplicación de ambos 
modelos, los resultados son muy similares. Es 
importante aclarar que la comparación se realiza 
para un suelo franco. Conforme la textura del 
suelo sea más arcillosa, es de esperarse que la 
diferencia en los tiempos de cómputo sea aún 
mayor. 

Otra ventaja que tiene el modelo simplifica- 
do sobre el modelo hidrodinámico completo es 
que la información necesaria para su aplicación 
es, en general, más fácil y menos costosa de ob- 
tener que aquella que se ocupa en la aplicación 
del modelo hidrodinámico completo. 


Diseño del riego por melgas 


Para el diseño del riego por melgas se pueden 
utilizar dos métodos (Rendón et al., 1997): 
pruebas de riego y modelos físico-matemáticos. 
Las pruebas de riego consisten en hacer expe- 
rimentos de campo, en donde las variables 
experimentales son el gasto de riego, longitud 
de la melga y tipo de suelo. Con base en los 
resultados de los experimentos, para cada tipo 
de suelo se obtienen combinaciones de gasto de 
riego unitario y longitud de la melga que hacen 
posible aplicar con un alto grado de uniformi- 
dad una lámina de riego específica. Las pruebas 
de riego tienen la desventaja de requerir un 
tiempo considerable para su ejecución, son cos- 
tosas y debe realizarse una cantidad elevada de 
pruebas para determinar el gasto de riego para 
el que se obtiene la mayor uniformidad en el 
riego; por estas razones son recomendadas sólo 
para calibrar los métodos de diseño. 

En riego es indispensable distinguir al 
menos tres eficiencias relacionadas con el pro- 
ceso: la eficiencia de aplicación, la eficiencia de 
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Figura 2. Frentes de avance obtenidos mediante la aplicación del modelo hidrodinámico completo (línea continua) y el modelo 


simplificado (línea punteada). Suelo franco de Montecillo. 
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Figura 3. Evolución en el tiempo del tirante y la lámina infiltrada en la cabecera de la melga, obtenidos mediante la aplicación 


del modelo hidrodinámico completo (línea continua) y el modelo simplificado (línea punteada). Suelo franco de Montecillo. 


requerimiento de riego y la eficiencia de uni- 
formidad del riego. La eficiencia de aplicación 
(Ea) se define como Ea = Vr/Vp, donde Vr es 
el volumen requerido para satisfacer las nece- 
sidades de agua en la zona de raíces del cultivo 
[L*] y Vp es el volumen de proyecto [L*]; Vr se 
obtiene con la expresión Vr = LnAr, donde Ln es 
la lámina de riego neta [L], definida conforme 
a los requerimientos de riego del cultivo, y Ar 
es la superficie de riego considerada [L?]. El 
volumen de proyecto se calcula como Vp = QpTr 
=LbAr, donde Qp es el gasto de proyecto [L*T*], 
Tr es el tiempo de riego [T] y Lb es la lámina de 
riego bruta. Si en la expresión (18) el numerador 
y denominador se dividen entre Ar se obtiene 
En = Ln/Lb. La eficiencia de requerimiento de 
riego (Er) se define como Er = Vd/Vr, donde Vd 
es el volumen disponible [L*]. Esta eficiencia 
indica la manera en que se están satisfaciendo 
las necesidades de agua del cultivo. 

La situación ideal respecto a la uniformidad 
se presenta cuando todas las plantas reciben la 
misma cantidad de agua, situación que equivale 
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a aplicar una lámina de riego uniforme en toda 
la longitud de la melga. Para evaluar la unifor- 
midad en la distribución de la lámina infiltrada 
puede utilizarse el coeficiente de uniformidad 
de Christiansen (CUC): 


CUC=1-+E (12) 


donde I, es la lámina infiltrada en el punto 1 
[L]; Í es la lámina infiltrada media [L], y nesel 
número de puntos considerados para efectuar 
el cálculo. 

Empleando un modelo numérico desarro- 
llado para el acoplamiento de las ecuaciones 
de Saint-Venant, y Green y Ampt, es posible 
determinar el gasto para el cual se obtiene la 
mayor eficiencia de uniformidad, manteniendo 
los valores lo más elevado posible en cuanto a 
las eficiencias de aplicación y requerimiento de 
riego; es decir, para determinar el gasto óptimo 
de riego. 
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Para calcular el gasto óptimo para diversos 
tipos de suelo, diferenciados según su caracte- 
rización hidrodinámica, se ha procedido como 
a continuación se indica. El valor de la presión 
en el frente de humedecimiento M), así como 
la conductividad hidráulica a saturación (K_) se 
determinan con base en los valores proporcio- 
nados por Rawls y Brakensiek (1981) en relación 
con la textura del suelo. El contenido volumé- 
trico de agua a saturación ha sido asimilado a 
la porosidad total del suelo (p), la cual ha sido 
determinada mediante el mismo procedimiento. 

El contenido volumétrico de agua necesario 
para asignar la condición inicial O =0, para la 
ecuación de Green y Ampt se determinó tenien- 
do en cuenta la humedad aprovechable de cada 
tipo de suelo, suponiendo que se ha consumido 
50% de dicha humedad antes de la aplicación 
del riego. La humedad aprovechable se deter- 
minó restando los contenidos volumétricos de 
agua correspondiente a la capacidad de campo 
(CC) y al punto de marchitamiento permanente 
(PMP), los valores de CC y PMP para cada tipo 
de suelo se estimaron de acuerdo con los repor- 
tados por Rawls y Brakensiek (1981), asociados 
con el triángulo de texturas. 

En el cuadro 1 se muestran los valores de los 
parámetros de la ecuación de Green y Ampt pa- 
ra diez tipos de suelo. El conocimiento de dichos 
valores y del valor del contenido volumétrico 
de agua inicial hace posible calcular la relación 


entre la longitud de la melga y el gasto de riego 
óptimo para cada suelo, dada una lámina de 
riego. 

El coeficiente de uniformidad varía de forma 
sensible con el gasto de riego. En la figura 4 se 
muestra esta situación para el caso de tres tipos 
de suelo diferentes. Los puntos de inflexión de 
cada una de las curvas coinciden con el valor del 
gasto óptimo para una longitud de melga de 100 
m y una lámina de riego de 10 cm para cada uno 
de los tipos de suelo especificados en la gráfica. 

Es posible determinar el gasto de riego óp- 
timo para diferentes longitudes de melga. La 
representación gráfica de dicho cálculo para un 
suelo franco se muestra en la figura 5 para tres 
láminas de riego diferentes. Puede apreciarse 
que existe una relación básicamente lineal entre 
el gasto óptimo de riego y la longitud de la mel- 
ga. El valor de la constante de proporcionalidad 
entre dichas variables puede expresarse en tér- 
minos de gasto de aporte por unidad de área, 
es decir, por unidad de ancho y por unidad de 
longitud de melga; el resultado se muestra en 
el cuadro 2 para tres láminas de riego típicas de 
8, 10 y 12 cm. El cuadro presentado puede utili- 
zarse con fines de diseño del riego por melgas. 


Conclusiones 


Se desarrolló un método de diseño de riego 
por melgas, obtenido mediante la aplicación de 


Cuadro 1. Parámetros medios de la ecuación de infiltración de Green y Ampt (1911) utilizados para el diseño del riego por 


gravedad. 
Suelo 9, (cm*/cm?) 8, (cm*/cm?) h, (cm) K. (cm/h) 
Franco arenoso 0.16 0.46 12 2.9 
Franco 0.20 0.46 25 LS 
Franco limoso 0.17 0.55 30 1.0 
Limo 0.14 0.50 50) 0.8 
Franco arcillo arenoso 0.18 0.42 12 2.0 
Franco arcilloso 0.25 0.48 38 0.4 
Franco arcillo limoso 0.26 0.49 60 0.15 
Arcilla arenosa 025 0.42 23) 0í5 
Arcilla limosa 0.32 0.48 100 0.05 
Arcilla 0.36 0.49 100 0.05 
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un modelo basado en el uso de las ecuaciones 
de Saint-Venant para el flujo del agua sobre el 
suelo, y la ecuación de Green y Ampt para el 
proceso de infiltración del agua en el suelo. La 
principal ventaja del modelo presentado es el 
poco tiempo de cómputo que requiere para su 
aplicación, comparado con un modelo hidro- 
dinámico completo reportado en la literatura, 
el cual emplea las ecuaciones de Saint-Venant 
para el flujo del agua sobre el suelo, acopladas 
internamente con la ecuación de Richards, que 
permite modelar el flujo del agua en el suelo. 
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Figura 4. Variación de la eficiencia de uniformidad con el gasto de riego para tres tipos de suelo. Lámina de riego de 10 cm. 
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Cuadro 2. Cuadro de diseño del riego por melgas: gasto en 1/s/m? para la aplicación Óptima de la lámina de riego. Factor de 
fricción k = 1/54; pendiente topográfica J,= 0.002. Modelo de ecuaciones de Saint-Venant, y Green y Ampt. 


Lámina de riego de 8 cm Lámina de riego de 10 cm Lámina de riego de 12 cm 
Suelo Qopt TR CUC Qopt TR CUC Qopt TR CUC 
1/s/m? horas - 1/s/m? horas - 1/s/m? horas - 
Franco arenoso 0.01565 1.6 0.921 0.01425 2.2 0.938 0.01330 2.8 0.951 
Franco 0.00965 2.6 | 0.922 0.00875 39 0.939 0.00810 a | ME 
Franco limoso 0.01040 2.4 0.922 0.00890 3.4 0.936 0.00800 4.5 0.949 
Limo 0.00810 3.1 0.924 0.00700 4.3 0.937 0.00620 5.9 | 0.950 
Franco arcillo arenoso 0.01740 1.5 0.923 0.01515 2.0 0.938 0.01360 2.7 0.951 
Franco arcilloso 0.00330 75| 0.927 0.00285 10.6 0.941 0.00270 13.4 | 0.952 
Franco arcillo limoso 0.00175 13.7 0.924 0.00145 20.9 0.940 0.00130 27.9 0.953 
Arcilla arenosa 0.00290 8.5 0.929 0.00265 11.4 0.945 0.00245 14.7 0.959 
Arcilla limosa 0.00060 40.6 0.888 0.00051 59.0 0.901 0.00046 77.5 0.911 
Arcilla 0.00053 46.5 | 0.888 0.00045 67.3 0.901 0.00040 90.3 0.911 
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Figura 5. Relación entre la longitud de melga y el gasto de aporte óptimo para un suelo franco para tres láminas de 
aplicación: 8, 10 y 12 cm; K en cm/s. 
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